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> ■ ABSTRACT 

o ■ 

We model the infrared (IR) emission from the ring-Uke dust disk around the main- 
sequence (MS) star e Eridani, a young analog to our solar system, in terms of a porous 
^ dust model previously developed for the extended wedge-shaped disk around the MS 

Q>^ ' star (3 Pictoris and the sharply truncated ring- like disks around the Herbig Ae/Be stars 

^ ■ HR4796A and HD 141569A. It is shown that the porous dust model with a porosity of 

~ 90% is also successful in reproducing the IR to submillimeter dust emission spectral 
energy distribution as well as the 850 iim flux radial profile of the dust ring around the 
more evolved MS star e Eri. Predictions are made for future SIRTF observations which 
may allow a direct test of the porous dust model. 

O ' Subject headings: circumstellar matter — dust, extinction — infrared: stars — plane- 

» ' tary systems: protoplanetary disks — stars: individual (e Eridani) 
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1. Introduction 

Together with Vega, Fomalhaut {a PsA), and f) Pictoris, the nearby (distance to the Earth 
d ~ 3.2 pc) low-mass (m* ~ O.Sm©) K2V main-sequence (MS) star e Eridani constitutes the "Big- 
Four" family of prototypical "Vega-excess" stars, which are MS stars with infrared (IR) radiation in 
excess of what is expected from their stellar photospheres (Gillett 1986). The IR excess associated 
with "Vega-excess" stars is generally attributed to the thermal emission of dust grains orbiting the 
central star in the form of a disk or ring and heated by its stellar radiation. The dust is thought to 
be "second generational" in nature rather than "primordial"; i.e., the dust is considered as debris 
generated by collisions of larger objects such as cometary and/or asteroidal bodies rather than 
remnants left over from the star formation process (Backman &; Paresce 1993; Lagrange, Backman, 
&; Artymowicz 2000; Zuckerman 2001). 

In recent years, e Eri has received much attention since it is thought to be analogous to the 
solar system at an age about 10-20% of that of the Sun: (1) in comparison with the other three 
"Big-Four" prototypes which all are much more luminous A-type stars, e Eri is closer in type, mass. 
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and luminosity (L* fa O.33L0; Soderblom & Dappen 1989) to the Sun (see Greaves & Holland 2000); 

(2) e Eri has a ring-like structure of ~ 55 AU in radius (which is very similar to the Kuiper Belt 
in size) and a central cavity inside ~ 30 AU (which is equivalent in scale to Neptune's orbit) as 
revealed by its 850 //m image obtained using the SCUBA {Submillimetre Common-User Bolometer 
Array) camera at the James Clerk Maxwell Telescope (JCMT; Greaves et al. 1998). (3) as the 
oldest (with an age of ^ 800 Myr; Henry et al. 1996) "Big-Four" Vega-type star, the formation of 
planets in its disk may have completed since planet formation is thought to be complete by about 
10 to 100 Myr (Lissauer 1993); (4) a roughly Jupiter-mass planet with a period of ~6.9yr appears 
to be present inside the central cleared region of the e Eri disk, as suggested by the radial velocity 
measurements spanning the years 1980.8-2000.0 (Hatzes et al. 2000). 

To better understand the origin and evolution of protoplanetary dust disks and the formation 
process of planetary systems, we need to know the physical, chemical, and dynamical properties 
of their constituent grains. In previous papers (Li &; Greenberg 1998, Li & Lunine 2003a,b), 
we have shown that a simple porous dust model is successful in reproducing the spectral energy 
distributions from the near-IR to millimeter wavelengths of the disks around the Herbig Ac/Bc stars 
HR4796A and HD 141569A and the Vega-type star /3 Pictoris, including the 9.7 ^m amorphous and 
the 11.3 /im crystalline silicate features or the 7.7 fim and 11.3 /im polycyclic aromatic hydrocarbon 
(PAH) features. In this Letter we will study the dust IR to submillimeter emission of the disk 
around the more evolved star e Eri. The objectives of this Letter are four-fold: (1) we wish to infer 
the physical and chemical properties of the dust in the e Eri disk and to infer the disk structure; 
(2) wc wish to gain insights into the evolution of dust mass, size, morphology, and mineralogy at 
different stellar ages; (3) we wish to know how widely applicable is the porous dust model; if it can 
be shown that the porous dust model is "universally" valid for disks at different evolutionary stages 
and with different geometrical structures, it will be a valuable guide for interpreting future data sets 
obtained by the upcoming Space Infrared Telescope Facility (SIRTF); (4) we wish to make mid-IR 
spectral, imaging and broadband photometry predictions for the e Eri disk; these predictions can 
be compared to future SIRTF observations in order to test further the validity of the porous dust 
model. 



2. The Porous Dust Model 

We have discussed in detail in Li &: Greenberg (1998) and Li & Lunine (2003a, b) that the 
dust in protoplanetary disks most likely consists of porous aggregates of protostellar materials 
of interstellar origin, although there is no consensus regarding the degree to which the original 
interstellar materials have been processed in protostellar nebulae. Unless there are IR signatures 
(e.g. the 11.3 /xm crystalline silicate feature) suggestive of the thermal history of the dust available, 
it is difficult to determine the composition of the dust as a priori. Therefore, in modeling the IR 
emission of dust disks we usually invoke two extreme dust types: cold- coagulation" dust made of 
coagulated but otherwise unaltered protostellar interstellar grains, and ^^hot-nebula" dust made of 
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grains that are highly-processed in protostellar nebulae where silicate dust is annealed and carbon 
dust is destroyed by oxidization (sec Li & Lunine 2003a, b for details). Since the e Eri disk is too 
cold to emit at the mineralogical "fingerprinting" mid-IR features (e.g. the 9.7 and 11.3 //m silicate 
bands and the 3.3, 6.2, 7.7, 8.6, and 11.3 /nm PAH bands), we will just adopt the "cold-coagulation" 
dust model. It has been shown in Li k, Lunine (2003a,b) that the major differences between the IR 
emission from these two dust models lie in their mid-IR spectral features, while their continuum 
emission are very similar. 

For the dust properties, we take the porosity (or fluffiness; the fractional volume of vacuum in 
a grain) to be P = 0.90 for 3 reasons: (i) dust models with P ~ 0.90 are successful in modeling 
the IR emission from the disks around HR4796A (Li & Lunine 2003a), HD 141569A (Li & Lunine 
2003b), and /3 Pictoris (Li k Greenberg 1998); (ii) a porosity in the range of 0.80 < P < 0.90 
is expected for dust aggregates formed through coagulation as demonstrated both theoretically 
(Cameron & Schneck 1965) and experimentally (Blum, Schrapler, & Kozasa 2003); (iii) a porosity 
of P ~ 0.90 for the "cold-coagulation" dust is consistent with the mean mass density of cometary 
nuclei for which the ice-coated "cold-coagulation" dust aggregates are plausible building blocks (see 
Greenberg & Li 1999). We assume a power-law dust size distribution dn{a)/da oc a~" which is 
characterized by a lower-cutoff flmirn upper-cutoff Omax and power-law index a; we take amin = 1 /^im 
and amax = 1 cm (see Li & Lunine 2003a,b).^ As far as the dust composition is concerned, we take 
the "cold-coagulation" dust model - the dust is assumed to be composed of amorphous silicate and 
carbonaceous materials (and H20-dominated ices in regions colder than ~ 110-120 K); the mixing 
mass ratios for the silicate, carbon and ice constituent grains are approximately derived in Li k 
Lunine (2003a), by assuming cosmic abundance, to be JTicarb/w-sii ^ 0.7 and fn\cc/ {fn^w + "icarb) ~ 
0.8 where msn, mcarb) and mice are respectively the total mass of the silicate, carbon, and ice 
subgrains. 

Guided by the 850 ;um SCUBA/JCMT image and the 850 ;um flux radial profile (see Figs. 1,2 
of Greaves et al. 1998), we will adopt a Gaussian functional formula for the dust spatial surface 
density distribution cr(r) = cxp exp[— 41n 2{(r — rp)/A}^] which describes the dust ring peaking at 
~ 55 AU from the star. The extent of the e Eri disk is cut off at an inner boundary of rin = 0.05 AU 
where small grains are heated to temperatures > 1500 K, and an outer boundary of rout = 200 AU 
outside of which there is little dust emission. 

We use Mie theory together with the Bruggman effective medium theory (Kriigel 2003) to 

calculate the absorption cross sections of the fluffy heterogeneous dust aggregates. Dielectric func- 
tions are taken from Draine k Lee (1984) for amorphous silicate dust; Li k Greenberg (1997) for 
carbonaceous dust; and Li k Greenberg (1998) for H20-dominated "dirty" ice. The e Eri stellar 
radiation is approximated by the Kurucz (1979) model atmosphere spectrum for K2V stars with 



^We assume all grains are spherical in shape; the grain size a is defined as the radius of the sphere encompassing 

the entire aggregate. The choice of amin and ttmax, as extensively discussed in Li & Lunine (2003a,b) for the disks 
around HR4796A and HD141569A, will not affect our conclusion. 
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an effective temperature of Teg = 5000 K, a surface gravity of lg5=4.0, and a solar metallicity 
Z/Zq = 1. The calculated dust IR emission spectrum is compared with the IRAS flux at 12, 25, 
60, and 100 /um and the SCUBA flux at 450 and 850 /^m (see Table 1 in Greaves et al. 1998). There 
also exist 1300 //m photometric measurements for the e Eri disk (e.g., Chini, Kriigel, & Kreysa 
1990, Chini et al. 1991), but we will not include these data since these measurements were all made 
using single-element bolometers with telescope beams smaller than the disk extent so that the flux 
obtained by these measurements may be underestimated (Zuckerman & Becklin 1993, Weintraub 
L Stein 1994). 

3. Model Results 

The model has four free parameters: a - the dust size distribution power-law exponent, - the 
disk radial location where the dust spatial distribution peaks; A - the FWHM of the dust spatial 
distribution; cTp - the mid-plane surface density at r = rp. Observational constraints include the 
observed IR spectral energy distribution and the SCUBA 850 /xm image. The latter is particularly 
important for constraining the spatial distribution of the dust in the e Eri disk. Instead of using the 
smoothed 850 /vm image presented in Greaves et al. (1998), we constructed an unsmoothed 850 /im 
image from all 1997-2000 data in the SCUBA archive using the data reduction procedure outlined 
in Bendo et al. (2003). However, both the image and the radial profile obtained from our analysis 
show an overall good agreement with those of Greaves et al. (1998). 

As shown in Figure 1, the P = 0.90 porous dust modeP with a ~ 3.1, rp ~ 55 AU, A ~ 30 AU, 
and fjp ~ 966 cm^^ provides an almost perfect fit to the observed IR to submillimeter emission 
(Fig. la) and to the azimuthally averaged flux radial distribution at 850 /xm (Fig. lb) except that 
the model is deficient in emitting at 850 /xm inside ~28AU (see below for further discussion).^ 

This model requires a total dust mass of ~ 4.22 x 10^^ g ~ 7.1 x 10^'^ m®. Let the goodness of 
fit be xV^^obs ^ Eilr {([Ai^A]^od - [Ai^Alobs) / [AAFA]„bs}' where [XF^U, and [XF^]^^^^ are 
respectively the observed and model-predicted flux densities, [AAF^J^j^^, is the observational flux 
uncertainty, A'obs = 6 is the number of data points to be fit. This model gives x^/^obs ^ 0.54. 

The apparent rise in flux density within r~4.5" (~15AU) was also seen in the smoothed 
SCUBA 850 //m radial profile of Greaves et al. (1998; see their Fig. 2). This rise may be associated 
with the small peak 5" -10" south of star'^ seen in the SCUBA 850 /um image (see Fig. 1 of Greaves 
et al. 1998). If this is real, we then need to resort to an additional dust component lying close to 
the star. This inner warm component would resemble the zodiacal dust cloud in our solar system. 



^In fact the porosity is reduced to P' 0.73 after ices fill in some of the vacuum (see Appendix B in Li & Lunine 

2003a) which occurs at r > 3AU. 

^The model-predicted 850 /xm radial profile is convolved with the SCUBA 14" beam using the actual SCUBA 
point spread function (PSF). 
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However, the SCUBA images are quite noisy in the inner region, so the anomalous emission inside 
r ~ 4.5" may be instrumental noise or image processing artifacts. We therefore exclude this "inner 
component" emission in our model. 

4. Discussion 

We have seen in §3 that the very same porous dust model shown to be successful in modeling 
the IR emission from the dust disks around the Herbig Ae/Be stars HR4796A (Li & Lunine 2003a) 
and HD 141569A (Li & Lunine 2003b) and the young MS star j3 Pictoris (Li & Greenberg 1998) 
is also successful in reproducing the IR to submillimeter SED as well as the 850 /xm flux density 
radial profile of the more evolved MS star e Eri (~800Myr old), although the model appears 
unable to provide sufficient 850 /um emission inside ^ 28 AU (see Fig. lb) . It is very likely that this 
inner component emission is not real (see §3); but if one really wants to push us to account for 
this emission, we find that it could be explained by invoking an inner "zodiacal" dust cloud of 
^ 5.2 X 10~^ mass located at ~ 2 AU from the star. 

SIRTF will be capable of sensitive imaging using the Infrared Array Camera (IRAC) at 3.6, 4.5, 
5.8, and 8.0 /xm, and using the Multiband Imaging Photometer (MIPS) at 24, 70, and 160 //m. SIRTF 
will also be able to perform low-resolution 5-40 /xm and high-resolution 10-37 /xm spectroscopic 
observations using the Infrared Spectrograph (IRS) instrument. However, the dust in the e Eri disk 
seems to be too cold to emit in an appreciable quantity in the IRAC bands (see Fig. la). Neither 
can this dust be heated sufficiently to emit at the characteristic mid-IR spectral features. But the 
MIPS imaging will provide powerful constraints on the e Eri dust spatial distribution. We have 
therefore calculated the MIPS 24, 70, and 160 /xm flux density radial distributions for the best-fit 
porous dust model (see Fig. Ic). All profiles are convolved with the MIPS model PSFs. In Table 1 
we show the model-predicted band-averaged intensities. 



Table 1: Dust IR emission (Jy) integrated over SIRTF bands predicted for the porous dust model. 
Also tabulated are the stellar photospheric emission (Jy). 



Instrument 


IRAC 


IRAC 


IRAC 


IRAC 


MIPS 


MIPS 


MIPS 


Aeflf 


3.6 /xm 


4.5 /xm 


5.8 /xm 


8 /xm 


24 /tm 


70 /tm 


160 /xm 


dust 


1.35 X 10-^ 


2.03 X 10-5 


3.43 X 10-5 


8.47 X 10-5 


0.10 


1.55 


1.71 


stellar 


65.4 


40.7 


28.6 


16.4 


2.26 


0.25 


0.056 



The dust in the e Eri disk is stable against the radiation pressure. As shown in Figure 2a, 
/3rp, the ratio of radiative pressure (RP) force to gravitational force is found to be smaller than 0.3 
for all grains. But the Poynting-Robertson (PR) drag would remove grains smaller than ~ 270 /xm 
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at r = 55 AU in a timescale (tpr) shorter than the age of the e Eri system (see Fig. 2b). ^ By 
integrating the PR dust removal rate (Att/S) pa^ /Tp^{a,r) [where p is the mass density of the dust 
aggregate: p ~ 0.25 (0.46) g cm~^ for the P = 0.90 {P' = 0.73 ice-coated) cold-coagulation type 
dust; see Appendix B in Li &: Lunine 2003a] over the whole size range and over the entire disk, we 
estimate the PR dust mass loss rate to be 4.3 x 10~^^ m^yr"^ for the porous dust model. Over 
the life span of e Eri, roughly 0.034 m® of dust is lost by the PR drag. Therefore, the observed 
grains need to be continuously replenished, very likely from cascade collisions and evaporation of 
larger bodies such as comets and/or asteroids. 

The inner "zodiacal" cloud, if it indeed exists, could be maintained by grains released from 
the sublimation of bombarding comets near the star.^ Taking the total dust mass loss ~3 x 10^^ g 
in the 1997 apparition of comet Halc-Bopp (Jcwitt & Matthews 1999) as a representative dust 
supply rate (per comet per apparition), we estimate a cometary infalling rate of ~ 7yr~^, roughly 
consistent with the estimation of <15yr~^ derived by Dent et al. (1995) from the observed CO 
upper limit assuming sublimating comets as a CO replenishment source in the e Eri disk.^ 

The grain-grain collision timescale is only 0.1 Myr for dust at Vp = 55 AU, significantly 
shorter than the PR timescale tpr. But it is unlikely for grain-grain collisions to effectively remove 
the dust from the disk; instead, their major role is to re-distribute the dust over different size bins 
through fragmentation. The partially cleared central cavity is most likely created by the accumu- 
lation of gas and grains into planctcsimals, as supported by the non-detection of circumstellar gas 
in this disk (Dent et al. 1995, Liseau 1999); and by the fact that the PR drag alone is not sufficient 
to clear the dust out to large radii up to ~30AU (Jura 1990, Greaves et al. 1998, Dent et al. 2000). 

Dent et al. (2000) modeled the e Eri SED in terms of compact grains, assuming the dust 
absorption efficiency to be Qabs = 1 for A<30//m and Qabs = (A/30/xm)"°-^ for A>30 fim, and 
adopting a uniform surface density distribution [i.e. a{r) = constant] in the range of 50<r<80 AU. 
Satisfactory fit to the observed SED was achievable. Using emissivities calculated from real dust 



Closer to the star, the threshold grain size below which grains will be ejected increases since tpr oc (see Eq.22 
in Li & Lunino 2003b). 

®In contrast, the cold outer ring is more likely to be coUisionally replenished, although the sublimation of cometary 
CO volatiles (with a sublimation temperature of ~20K) occurring inside r<110 AU does generate dust outflow. But 
dust outflow only becomes significant when comets are close to the star, say, at r < 3AU where H2O ice starts to 
sublimate. 

®If we adopt the Halley dust loss rate ~1 x 10^'' g per orbit (Lisse 2002), we would expect a comet infall rate 
~ 300 times more frequent. This is much higher than what is observed (~ 20-30 comets per year) for the present 
solar system. But we note that the number of cometary bodies and the bombardment rate were much higher at the 
early stage of our solar system. The solar system zodiacal cloud has a mass loss rate of ~ 5 x 10"^'' yr"^; Lisse 
(2002) argued that the short-period comets alone are enough to supply the dust. In the (3 Pictoris disk, the number 
of observed cometary bodies (crossing the lino of sight) is ~ 300-400 per year (Beust &: MorbidoUi 1996); however, the 
actual number would be much higher since many objects do not cross the line of sight and are therefore not visible 
in spectroscopy. 
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mixtures, however, Sheret, Dent, & Wyatt (2003) found that compact grains are not able to fit the 
observed SED. 

The model presented in this Letter is somewhat simplified in the sense that we have assumed a 
symmetrical morphology for the e Eri disk. Asymmetries and discrete flux enhancements ("clumps" 
or "bright spots") arc actually seen in its 850 /um SCUBA images (Greaves et al. 1998). As suggested 
by Ozernoy et al. (2000) and Quillen & Thorndike (2002), these features may be caused by the 
capture of grains into the mean motion resonances with a planet of Jupiter-mass embedded in the 
disk. This simplification would not affect the general conclusions of this Letter since the IR emission 
is mainly determined by the general distribution of dust in the disk, that is, a narrow ring. 

5. Conclusion 

The porous dust model, previously developed for the extended wedge-like disk around the 
~ 15 Myr-old MS star f3 Pictoris and the narrowly confined ring-like disks around the Herbig Ae/Be 
stars HR4796A (~8Myr old) and HD 141569A (~5Myr old), is shown also applicable to the dust 
disk around the ~ 800 Myr-old MS star e Eri, the oldest "Big- Four" Vega-excess star. Modeled 
as a ring peaking at ~ 55 AU from the star with a FWHM ~ 30 AU, analogous to the young solar 
system Kuiper Belt, the observed IR to submillimeter emission as well as the 850 /itm flux density 
radial distribution of the e Eri disk is closely reproduced by porous grains of ~ 7.1 x 10~^ mass 
consisting of --^90% vacuum in volume. Predictions are made for future SIRTF/MIPS imaging 
observations that will offer a useful test of the porous dust model. 
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for helpful discussions and/or suggestions. A. Li thanks the University of Arizona for the "Arizona 
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by a grant from the NASA origins research and analysis program. G.J. Bcndo is a Guest User of 
Canadian Astronomy Data Centre, which is operated by the Dominion Astrophysical Observatory 
for the National Research Council of Canada's Herzberg Institute of Astrophysics. 
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Fig. 1. — (a) IR to submillimeter emission of the e Eri system (open circles: the 12, 25, 60, and 100 /^m IRAS 
data and the 450 and 850 SCUBA data) compared with the model spectrum (solid line; which is the sum of 
the stellar [dashed line] and the dust [dotted line] contributions), (b) Radial profile of the dust emission at 850 /im. 
Open squares plot the azimuthally averaged flux density radial distribution calculated from the SCUBA 850 fim 
image (also see Figs. 1,2 of Creaves et al. 1998 but theirs were smoothed with an 8" Caussian). The apparent rise 
within r ~ 15 AU is very likely artificial (see §3). Solid line plots the model- predicted 850 /im emission radial profile 
convolved with the SCUBA 14" beam, (c) Radial profiles predicted for the SIRTF MIPS bands at 24 /^m (solid line), 
70 fj,m (dotted line), and 160 /xm (dashed line). All profiles are convolved with the MIPS model PSFs and normalized 
to their peak values: « 0.71, 6.47, and 6.12 mjy arcsec"^ for the 24, 70, and 160 /im bands, respectively. 
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Fig. 2. — (a) Ratio of the radiative repulsion to the gravitational attraction (/3rp) for the best-fit "cold-coagulation" 
dust (without ice P — 0.90 [solid line] or with ice P' = 0.73 [dashed line]), (b) The orbit decay timescales tpr due 
to the Poynting-Robertson drag for the best-fit "cold-coagulation" dust (without ice P — 0.90 [solid line] or with ice 
P' = 0.73 [dashed line]) at a radial distance of rp = 55 AU from the central star (note rpR, oc r^: the PR timescale 
decreases for dust at a smaller radial distance). The dot-dashed horizontal line plots the e Eri age (« 800 Myr). 



